


Institutional Archive of the Naval Postgraduate School 


Calhoun: The NPS Institutional Archive 
DSpace Repository 


Theses and Dissertations 1. Thesis and Dissertation Collection, all items 


1973 


An optimum low frequency digital 
communications receiver 


Doherty, James Thomas 


: Massachusetts Institute of Technology 


http://ndl.handle.net/10945/16733 


Downloaded from NPS Archive: Calhoun 


| Calhoun is the Naval Postgraduate School's public access digital repository for 
D U DLEY research materials and institutional publications created by the NPS community. 
sa Calhoun is named for Professor of Mathematics Guy K. Calhoun, NPS'‘s first 
KNOX appointed — and published — scholarly author. 


LIBRARY Dudley Knox Library / Naval Postgraduate School 
411 Dyer Road / 1 University Circle 
Monterey, California USA 93943 





http://www.nps.edu/library 


AN OPTIMUM LOW FREQUENCY DIGITAL 
COMMUNICATIONS RECEIVER, 


James Thomas Doherty 















‘POSTGRADUATE SCHOBD 
2EY, CALIF. 93940 


ae 
AN OPTIMUM LOW FREQUENCY DIGITAL COMMUNICATIONS RECEIVER 


by 
James Themes Doherty Jr. 
Y 


B.S., U.S. Coast Guard Academy 
(1969) 


SUBMITTED IN PARTIAL FULFILLMENT OF THE 
REQUIREMENTS FOR THE DEGREES OF 
ELECTRICAL ENGINEER 
AND 
MASTER OF SCIENCE 


at the 


MASSACHUSETTS INSTITUTE OF TECHNOLOGY 
June 1973 





LIBRARY 
NAVAL POSTGRADUATE SCHOOD 
MONTEREY, Ca: Tr. 93940 


=e 
AN OPTIMUM LOW FREQUENCY DIGITAL COMMUNICATIONS RECEIVER 


by 


James Thomas Doherty Jr. 


Submitted to the Department of Electrical Imgineering 
on 11 May 1973 in partial fulfillment of the require- 
ments for the Degrees of Electriztal Ingineer and Master 
of Science, 


ABSTRACT 


In the Low Frequency (LF, 30-300 kHz) band of the 
electromagnetic spectrum, one finds many digital commu- 
nication systems. At these frequencies, atmospheric 
noise, predominantly caused by radiation from lightning 
discharges, is the dominant noise which corrupts the sige 
nal. When viewed at the output of the receiver's band- 
limiting filter, the noise waveform is very definitely 
non-saussian. Since the noise is non-gaussian, it would 
be difficult to design an optimum communication receiver 
without knowledge of specific properties of the noise. 

A model for bandlimited low frequency atmospheric 
noise, which accurately fits both the first order 
probability density and the time structure of observed 
noise waveforms, and which is well suited for digital 1 
computer simulations, has been developed by Feldman. 

In the development of this model, a statistical dependence 
between the noise waveform on one bandlimited freauency 
channel and the envelope of the noise on a disjoint, but 
nearby, frequency channel was identified and incorporated 
into the model. On the basis of this property, a design 
for a two channel low frequency communication receiver, 
which makes nearly optimum use of the joint channel 
information, was proposed. 


THESIS SUPERVISOR: Amar G. Bose 
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This work has addressed itself to the task of 
Simulating the joint channel receiver and testing its 
performance relative to a single channel linear corre- 
lation receiver. Both receivers operate in the presence 
of low frequency atmospheric noise, as simulated by the 
above mentioned noise model. 

The modem was assumed to be coherent phase shift 
keying, with the signal set consisting of + A cos oni o i; 
The criterion for performance comparison was an estimate 
of probability of error per binary bit of information 
received. Our results show that, with the dual channel 
receiver, performance improvements of more than an order 
of magnitude can be expected during periods of severe 
non-gaussian noise activity. Svegestions for further 
research with this receiver design are presented in the 
last chapter of this work. 
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Chapter 1 
INTRODUCTION 


In the Low Frequency (LF, 30-300 kHz) band of the 
electromagnetic spectrum, one finds many digital commu- 
nication systems. At these frequencies, the dominant 
noise encountered is atmospheric noise, which is 
predominantly caused by radiation from lightning 
discharges. Since individual discharges tend to 
dominate the noise waveform observed at the output 
of a bandlimiting filter, we must assume that the noise 
cannot be statistically described as a gaussien noise 
process. Clearly then, one cannot hope to build an 
optimum communication receiver without some knowledge 
of the specific non-gaussian characteristics of the low 
frequency atmospheric noise. 

The research reported in this work will demonstrate 
one particular realization of a near optimum binary 
digital communication receiver which is based upon the 
observed characteristics of bandlimited atmospheric 
noise. In this chapter, we will discuss the actual 
atmospheric noise and a model for that noise which 
can be used in digital computer simulaticns, 

In Chapter 2, we present a receiver design which 
uses a second disjoint frequency channel to provide more 


information about the noise which corrupts the signal. 
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We will show that this receiver makes nearly optimum use 
of the joint channel information. We will also compare 
this receiver to a single channel linear correlation 
receiver. 

Chapter 3 is devoted to an explanation of the 
calculations necessary to relate the analog system to 
the digital simulations of the proposed receiver. Our 
criterion for evaluation of receiver performance will 
simply be an estimate of the probability a peneer per 
bit of binary information received. 

In Chapter 4, we will present the results of our 
Simulations of both the single channel and joint channel 
receiver models, in the presence of the simulated low 
frequency atmospheric noise. We will give bounds on 
the accuracy of our results and explain some of the 
problems encountered in carrying out the simulations. 

Finally, Chapter 5 will give a summary of this 
research, the conclusions reached, and some suggestions 


for further research. 


1.1 Low Frequency Atmospheric Noise 


There have been many observations of bandlimited 
low frequency atmosvheric noise, which have led to 
descriptions of the noise characteristics (for example, 
see 1,2), but all note that the noise waveform is an 


extremely non-gaussian process. At individual times, 
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the noise waveform is dominated by peaks of high 
instantaneous energy, with several peaks occasionally 
occurring rapidly enough to yield a seemingly continuous 
burst of high energy. 

Fach high energy peak can be related to the 
radiation released from a lightning discharge. There 
are two basic strokes in each discharge: the leader 
stroke with an average length of about 1 msec and 
centered in frequency about 30 kHz, and the return 
stroke with a length of about 100-200 usec and centered 
about 3-8 kHz. In many instances, there are complicated 
repititions of these two basic strokes which yield the 
continuous noise burst behavior noted above. Between 
high energy bursts, the noise waveform can be accurately 
characterized as a low level a0 este waveform, 

Pelaneee reported that the following observed 
properties of the bandlimited noise waveform were useful 
in developing a model for this noise: First, at low 
frequencies, the noise is highly dependent upon local 
weather conditions. Second, above a certain threshold, 
the first order joint statistics of the noise waveform 
in one frequency channel and the noise envelope in a 
second disjoint, but nearby, frequency channel are 
Statistically dependent. Third, the bandlimited waveform 
has an average power spectral density which is proportional 


to the square of the narrow-band filter's frequency 
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response. And fourth, for bandwidths of 1 kHz or larger, 
the process can be considered uncorrelated beyond 1 msec. 
In the next section, we will present the noise model which 
Feldman subsequently developed. 

In his analysis of atmospheric noise, enemas noted 
that one of four basic weather/noise conditions is 
sufficient to describe the noise at any particular time. 
These noise conditions are quiet, tropical, frontal, and 
quiet night. Quiet corresponds to a stable air mass, 
with little thunderstorm activity, in the vicinity of the 
receiver. Tropical conditions exist when warm, moist, 
unstable air, with much thunderstorm activity, is present 
near the receiver. Frontal corresponds to the passage 
of a front in the vicinity of the receiver. Quiet 
night represents a transitional period between quiet and 
tropical conditions. The order of severity varies from 
nearly gaussian to nearly continuous high energy burst 
behavior for quiet, quiet night, tropical, and frontal 


conditions, respectively. 


1.2 Noise Model 
The model for bandlimited low frequency atmospheric 
noise, y(t}, which we will consider in this work is the 
Por Vovaee 
y(t) = ny(t) + a(x(t))no(t), (1.2) 


where n,(t) and no(t) are statistically independent 
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gseussian processes and ny (t) represents the low level 
gaussian background process. The quantity a(x(t))no(t) 
represents the high energy bursts of noise due to 
lightning discharges, and x(t) is a two-state Markov 
process with a time-varying rate parameter. The two 
states of x(t) turn on and off the high energy noise 
component. Whenever x(t) is in the "off" state, a(.) 
is zero, and whenever x(t) is in the "on" state, a(.) 
becomes a random number, which is constant over the 
time interval that x(t) remains in the "on" state. 

Each time x(t) changes from "off" to "on," a new random 
value for a(.) is generated, and a(.) remains constant 
until x(t) turns it off. 

This model incorporates the observed dependence 
between the noise waveform on one channel and the 
envelope of the noise on a second disjoint, but nearby, 
frequency channel by assuming that a(x(t)) remains 
independent of frequency. Thus, it follows that the 
model for the envelope of the bandlimited noise, v(t), 
on a second frequency channel is 

v(t) = vy (t) + a(x(t))vo(t), (eZ } 
where v,(t) and vo(t) are statistically independent 
Rayleigh processes. 

It must be noted that this noise model cannot 


exactly specify the short-time structure of the actual 
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atmospheric noise. However, since LF receivers use sharp 
bandlimiting filters, with long observation times, this 
noise model is sufficient because it accurately simulates 
the average behavior of the noise over these relatively 
long tine intervals. The model is based on the first 
order probability density of the bandlimited noise. In 
addition, properties of a(x(t)) are such that the model 
accurately simulates the actual long-time time structure 
of the noise. 

Each of the four noise conditions (quiet, quiet night, 
tropical, and frontal) can be simulated with the noise 
model simply by varying the parameters used in the 
generating mechanism. In this research, all four noise 
conditions were simulated extensively using a canonic 
noise generator and table of associated parameters 
Siven in Figure 3-11 and Table 3-1, respectively, of 
reference (1). Since these two items are of vital 
importance to this research, they have been included 
as Figure 1.1 and Table 1.1 of this report. The results 
of the simulations, conducted in this research project, 


are given in Chapter 4 of this work. 
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Figure 1.1 Canonic Noise Generator 
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Table 1.1 Noise Model Parameters 


(braces and arrows indicate parameters we used) 





Chapter 2 


RECEIVER 


In Chapter 1, we presented some properties of the 
bandlimited low frequency atmospheric noise, the noise 
which corrupts signals in low frequency redio communication 
systems. We also presented a model, suitable for simulation 
on a digital computer, which accurately represents that 
noise. The low frequency atmospheric noise, when viewed 
at the output of the receiver's bandlimiting filter, has 
been observed to be extremely non-gaussian. However, we 
presented several properties of the noise which can be 
exploited in the design of a nearly optimum receiver. 

In this chapter, we will present a model for a 
phase coherent linear correlation receiver. We will 
show that this is the optimum receiver for signals 
corrupted by white gaussian noise, but it is not neces- 
Sarily optimum when the noise is non-gaussian. We will 
also present a receiver design, based on the specific 
observed properties of bandlimited low frequency 
atmospheric noise, which should yield nearly optimum 
performance. However, implementation of the new design 
ls approximately twice as complicated as implementation 


of the linear correlation receiver, 
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2.1 Modulation-Demodulation Schene 

Since the basic purpose of this research is to 
compare the performance of a new receiver design with 
that of a standard receiver design, we have chosen a 
computationally simple modulation-demodulation scheme. 
We postulate, since we are comparing both receivers 
under the same set of conditions, that our results in 
this somewhat idealistic case will still be representative 
of the results we would expect in a more realistic case. 

In the simplest case, the receiver must decide 
whether or not a signal is present. In a slightly more 
complicated, and more realistic, case, the receiver must 
decide between two signals. We have thus chosen our sig- 
nal set to contain two signals, each containing the same 
amount of energy. This signalling scheme, in which one 
Signel is the negative of the other, is known as anti- 
podal Seen 

Since only a limited amount of data was analized in 
the development of the noise model, we have parameters 
avallable to model the bandlimited noise waveform only 
under certain frequency and bandwidth constraints. 
Keeping these constraints in mind, we must satisfy the 
not unrealistic condition that all radiated signal 
energy be confined to the frequency band centered at 


65 kHz, with a bandwidth of 1 kHz.” (see Table 1.1) 
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This means that in the baseband, the lowpass region of 
the frequency spectrum, our signals must contain no 
energy in any frequency component greater than 1 kHz. 
In order to satisfy this constraint in a simple menner, 
we have chosen our baseband signal set to consist of 

Sp(t) = + Ay (24 ) 
and thus, since A is a constant, the frequency spectrum 
of each of our baseband signals consists only of an 
impulse at the frequency origin. 

However, we are restricted to a frequency band 
centered not at the origin but at 65 kHz for transmission, 
Tnerefore, we must modulate our signals up to this pass- 
band. To do this, we simply allow our signals to ride 
on a cosine carrier of frequency 65 kHz; and thus our 
passband signal set becomes 

s(t) = + A cos nf ats (22 ) 
where f, = 65 kHz. 

At this point, we will separate the members of our 
transmitted signal set, thus: 

So(t) = A cos 2nf it, (2.3a) 
which represents a binary zero being sent, and 
s,(t) = -A cos 27f,t, (2.3b) 
which represents a binary one being sent. From the above, 
it is clear that 
S,(t) = “So(t) = A cos (27mf,t + 8), (2. 3c) 





— 


where @ =7Tl. Thus, we see that our signal set consists 
of two phase shift keyed (PSK) signals. 

In general, this is not a practical way of signalling 
by radio communication systems due to the high voltage 
switching transients which occur when reversing phase. 
However, for signals corrupted by gaussian noise, the 
optimum PSK demodulator performs 3 db better than the 
optimum frequency shift keyed (FSK) demodulator, where 
binary FSK signal sets consist of 

Sq (t) = A cos 27fot (2.4a) 

and 

s,(t) = A cos 2nf,t, fy # fo. (2.4b) 
Furthermore, low frequency signalling schemes which 
approach PSK in performance have been developed, and it 
is felt that analysis of a PSK model will give an 
indication of performance to be expected in these new 
Peienes. 

If there were no noise present, the optimum demodu-e 
lator would filter out all signals but those in the 
passband, monitor the received signal, compare it with 
a reference signal (which is an exact replica of either 
Sy (t) or s,(t)), and decide which signal was sent. 
However, there invariably is noise present, which in the 
Simplest form is additive white gaussian noise. Other 
forms of noise, such as a random attenuation of the 


Signal or a random phase shift, may exist. We will, 
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however, restrict ourselves to tne additive noise case. 
The reasons for this will be given in the next section, 


where we explain the linear correlation receiver, 


2.2 Linear Correlation Recelver 

In almost any text on communication theory, one 
can find a description of a linear correlation receiver 
(see, for example, 3,4,7). In Figure 2.1, we show a 
linear correlation receiver which decides which of two 
arbitrary signals has been received. In Figure 2.2, we 
show how that receiver can be simplified for the case of 
antipodal signals, the case we are considering. 

Let us now assume that the received signal, r(t), 
consists of the following (see Figure 2.2): 

r(t) 


s(t) + n(t) 

= + A cos 2nf,t + n(t), (2.5) 
where n(t) is the additive white gaussian noise. We 
will assume that n(t) is a stationary, zero mean pro- 
cess, with autocorrelation function, 


RA (ty-to) = $(t,-to), (2.6) 


where §(.) 1s the Dirac delta function. From this, it 
ls clear that 


r,(t) = + AS cos“ (27f,t) + An(t) cos 21f,t 
= + AS (1 + cos 4i7f,t) + An(t) cos 27f,t, (2.7) 
2 


and 
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Figure 2.1 General Linear Correlation Receiver 
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Figure 2.2 Linear Correlation Receiver-Antipodal Signals 
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where we have assumed that T, the bit length or length 
of time in which the receiver must decide which signal 
was sent, is equal to an integral number of periods of 
the cosine wave. Actuelly, this is not a restrictive 
assumption because the bit length is typically much 
longer than the leneth of one period of the cosine wave. 

It is clear, since n(t) has been defined to be a 
gaussian random process, and since ro is the result of 
a linear filtering operation on n(t), that no iste 
gaussian randon variable.? Thus, we can completely 
specify the probability density of ro if we can find 
the first and second moments of i250 


Taking the expected value of equation 2.8, we get 


2 T 
E(ro) =+ AP +1 E( f An(t) cos 2nf,t dt) 
Zz fi: ‘OQ 
2 of: 
=e oe AEMtec) yecosment cd t 
Z ee 4 Cc 
= + A‘ ‘ (2.9) 
2 


Taking the second moment, we obtain 


T 
E(r2) = (+ a&)* + a2 af 2 An(t) cos 2nf,t dt) 
2 27 0 


eee 
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Substituting equation 2.6 into equation 2.10, using the 
properties of the Dirac delta function and trignometric 


identities, we obtain 


2 2,2 2 
2 5 aoe 
Therefore, 
Var(r,) = A“Noe- (2.12) 
2° yr 


Let us now define the following hypotheses: 


Ho? r(t) = so(t) + n(t), (2.13a) 
Hy: r(t) = sy(t) + n(t), (2.130) 
which can be reduced to 
a 2 
Hot Yo is N(A’, A“No) (2.14a) 
ae a he 
H,: Yo is N(-A%, A2N,) (2.14) 
fe aie. 


where N(m,n) denotes a gaussian random variable with menn 
mand variance n. In the case where either hypothesis 


ls equally likely, implying that either a binary one or 
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@ binary zero has been sent with equal probability, 
end where the penalty for making an incorrect decision 
is independent of whichever hypothesis is true, the 
decision rule for minimum probability of error is 
alwayss (see for example 3,4) 
Choose Hp if ro > 0; (2.15a) 
Choose H; if ro a70% (2a) 
Thus, we have shown that Figure 2.2 is a model of 
the optimum receiver for our signal set when the noise 
is geussian. However, we have not yet explained why we 
considered only additive noise. The concept of 
modelling the noise as inflicting a random attenuation 
uvon the information part of the received signal does not 
change the optimum receiver structure for the signal set 
we have chosen. All operations remain the same, and the 
decision threshold remains at the origin; however, the 
probability of making an error will increase. Since 
the receiver structure remains the same, we feel that 
there is no necessity for considering attenuation. 
However, the random phase shift is not easy to 
disregard. The primary cause of random phase shift noise 
is the fact that there is no way for our receiver model 
to insure that its reference oscillator will remain 
phase locked with the transmitter's oscillator. When 


unknown phase errors exist, the receiver loses its 
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ability to distinguish between Sg({t) and s;(t). The 
receivers we have shown in Figures 2.1 and 2.2 are 
known as coherent phase receivers since they assume 
that their oscillators are perfectly locked to the 
transmitter oscillator, In reality, receivers are 
phase incoherent, and they are considerably more 
complicated than the simple coherent receivers we have 
showmm. (see for example 3,4,5,6) We feel, however, 
Since the object of this research is to perform a 

first test uvon the dual channel receiver design, which 
we present in the next section, that we should only test 
the simple coherent phase receivers. Moreover, we feel 
that, since both receivers we are comparing will be 
coherent phase types, the results should be somewhat 
typical of the results to be expected in the incoherent 
phase case, Finally, if our results indicate signif- 
icantly improved performance, then oe ee dif fereul t, 
and more realistic, incoherent phase test would be 


justified. 


2.3 Joint Channel Receiver 


In the last section, we presented an optimum 
receiver for use when the corrupting noise is gaussian. 
However, when the noise is non-gaussian, we have no 
guarantee that an optimum gaussian receiver will yield 


optimum performance. With knowledge of the svecific 





Se 
nonegaussian nature of the actual noise encountered in 
low frequency radio communication systems, it should be 
possible to design a receiver which performs better than 
one which was designed for gaussian noise. 

Receiver designs which have been tried, or proposed, 
for dealing with the non-gaussian high energy noise 
bursts generally involve some type of nonlinear deen 
These receivers use clippers, hole punchers, and other 
similar nonlinear devices to limit the asiinatine effect 
of the high energy noise bursts, 

Feldmen* proposed another approach to the receiver 
design problem. He noted that, above a certain threshold, 
the noise waveform on the signal channel exhibited a 
statistical dependence upon the value of the noise 
envelope on a nearby, disjoint bandlimited frequency 
channel, His design exploits this characteristic in a 
nearly optimum manner by applying a multiplicative 
weighting function, based upon the inverse of the 
envelope measurement, to the signal channel. As it 
applies to this work, this receiver design, Figure 5-1 
in reference (1), is shown here as Figure 2.3. It 
should be noted that this receiver is basically a 
modified linear correlation receiver. 


The parameter v, is critical to this design, and 


t 
it is given by equation 4.32 of reference (1), which is 
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Figure 2.3 Dual Channel Receiver 
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repeated here; 
ve" em 6? aye / Onn 2s (2.16) 
where the quantities on the rignt hend side are given in 
Table 1.1 of this work. Whenever the noise envelope on 
the second channel, v(t), exceeds the threshold value, 
Vz. there exists a statistical dependence between v(t) 
and the bandlimited noise waveforn, y(t); on the 
information channel. Thus, a weighting function, z(t), 
Ziven below, is applied to the information channel: 
Zit edie aca ve 
v, / v(t), v(t) > v,. (Zi s 
We make no attempt to analytically verify that this in- 
verse first order weighting function is indeed nearly 
optimum. Instead, we rely upon our results to verify 
this fact. The weighting function given in equation 2.17 
above was first chosen as intuitively vleasing,. 

We verified our choice of weighting function 
experimentally. First, we point out that the object of 
the dual channel receiver is to reduce the problem from 
one in which we have a constant signal power level and 
a time varying, and unknown, (non-gaussian) noise power 
level to a problem in which the noise power level is kept 
constant and the signal power level is time varying in 
a known manner. Keeping this objective in mind, we 


conducted simulations for all four weather/noise 
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conditions, using various weighting functions. We 
Simultaneously simulated an ad hoc receiver which was 
composed of a time varying signal, formed by applying 
the weighting function to only the information part of 
the received signal, and a gaussian noise process. We 
then chose the weighting function which yielded dual 
channel receiver performance closest to the performance 
of the ad hoc receiver. The final weighting function, 


z(t), thus determined is given by equation 2.17 above. 


2.4 Chapter Summary 


In this chapter we have presented the modem to be 
used in the receiver simulations of this work. In 
addition, we have presented an analysis of a linear 
correlation receiver which operates on binary antipodal 
Signals that have been corrupted by additive white 
gaussian noise. We have also presented an improved 
design for a low frequency digital receiver, which 
utilizes information about the specific non-gaussian 
nature of the low frequency atmospheric noise to yield 
nearly optimum performance, 

In the next chapter, we will explain some of the 
problems encountered in attempting a digital computer 
Simulation of these receivers. We will also give examples 
of the calculations necessary to relate the digital 


simulation to the analog system being evaluated, 
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Chapter 3 


SIMULATION CALCULATIONS 


In this chapter, we will present the information 
necessary to relate our digital models of the noise 
and receivers to the actual analog systems we simulated. 
In the first section, we present the calculations 
necessary to convert the passband receiver model to 
baseband. In the second section, we calculate the 
probability of error for the simulated gaussian noise 
case. We also explain how we checked our receiver and 
noise model implementations on the computer. In 
Section 3.3; we give the calculations which relate the 
sampled data noise power spectrum to the analog noise 


power spectrum, in the white gaussian noise case, 


3.1 Baseband Calculations 
The received Signals which we described in the 


last chapter were located, in frequency, in the receiver's 
passband. It is a well known fact that a signal must 
be sampled at least twice as fast as the most rapidly 
varying frequency component in the signal, in order to 
be able to recover the signal from the samples. 

The highest frequency component of our received 
signal, r(t), given by equation 3.1, is the carrier 


frequency plus the bandwidth. The received signal, with 





SO ilu 

y(t) representing the bandlimited atmospheric noise, is 
Gt) eamcas Zita bat VC (3,1) 

In our case, the carrier frequency is 65 kHz, and the 

bandwidth is 1 kHz; thus, the highest frequency comyvonent 

allowed in r(t) is 66 kHz. Therefore, to accurately 

recover r(t) from the samples, we must sample the signal 

at a sampling frequency of at least 132 kHz. 

‘ Bit lengths for low frequency digital communication 
receivers, the length of time in which the receiver must 
decide which of the two possible transmitted signals has 
been received, are usually about 20-50 MES Oe For 
simulation purposes, we have chosen the lower end of 
this scale, in order to keep each simulation run as 
short as possible. Sampling at 132 kHz yields a time 
between samples of about 0.0076 msec, which we will round 
up to 0,01 msec for convenience. To simulate one bit of 
information, or 20 msec, at this samvling rate, we would 
require 2000 samples. 

Our evaluation of receiver performance is based 
strictly upon an estimate of the probability of an 
error being made by the receiver in any particular bit 
of binary information received. In order to form this 
estimate, we have chosen the following approach: count 
the number of bits simulated, count the number of incorrect 
decisions made by the receiver, and divide the latter by 


the former. Olsen” recommends that if the maximum number 
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of bits we simulate is 1000, then the minimum probability 
of error we can estimate, with any accuracy, is 0.01. 
Further details about, and statistical bounds upon, our 
estimation procedure will be given in Chapter 4, 

Based on the above, we feel that 1000 bits is the 
minimum number of bits we must simulate in order to 
derive meaningful results. In this case, to simulate 
the passband system, we will need 2,000,000 samples of 
the analog waveform, taken at a sampling interval of 
0.01 msec. This is incredible, especially when it is 
realized that 1000 bits represent only 20 seconds of 
Slenalling. This is also impractical for various 
reasons: the first reason is the tremendous computation 
time needed, and the second is that our noise model is 
incapable of accurately simulating such fine time 
structure of the noise process as is required by this 
sampling eee 

We have thus far ignored the one fact which allows 
us to solve this dilemma, the noise and transmitted 
Signal are bandlimited. Therefore, we can accurately 
convert the passband system to a baseband system in 
wnich the maximum frequency component is 1 kHz. Clearly 
then, the maximum sampling rate necessary is 2 kHz, and 
the time between samples is 0.5 msec. Therefore, the 
number of samples per bit (20 msec bit length) is only 


40, and the number of samples per 1000 bit simulation 
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run is only 40,000. Simulation of the baseband receiver 
is, of course, a much more tractable problen. 

Any bandlimited process, random or deterministic, 
can be reduced to lowpass (baseband) processes. The 
procedure for doing this is rather simple, and it is 
presented in various references on communication theory 
(see, for example 3,4). We simply multiply the passband 
process by a cosine and/or sine of the carrier (bandcenter) 
frequency, and pass the resulting process through a low 
pass filter. <A criterion for this procedure to work is 
that the bandwidth of the passband process must be much 
less than the carrier frequency. 

Since the procedure which converts a passband process 
to @ baseband process is actually a linear filtering 
operation, and since the received signal, r(t), is a linear 
combination of the transmitted signal, s(t), and the 
corruptive atmospheric noise, y(t), we can treat the 
conversion of each of these processes separately and 
linearly combine the results. 

First we one ider the transmitted signal (we consider 
eray s(t) = S(t), since s,(t) follows directly): 

s(t) = A cos 27f,t. ae 
Then the baseband process, by cosine conversion, s(t), 
before low pass filtering, is 


so(t) = A cos*(2nfgt) = A (1 + cos 4nf,t), (3.3) 
2 
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and, of course, after filtering, 
so(t) = A/2. (3.4) 
Similarly, by sine conversion, we get s(t): 


s(t) = A (cos 2nf,t)(sin 2rf¢t) 


‘ sin 4nf ts (3.5) 
and after low pass filtering, this becomes 

Seon eee (3.6) 
Therefore, we see that in order to convert the transmitted 
Ssignel portion of r(t) to a baseband process, we simply 
must calculate s,(t). 

We now consider the baseband conversion of y(t), the 
bandlimited low frequency atmospheric noise process. By 
cosine conversion, before low pass filtering, we obtain 
Y(t): 

nt) = (bt) fcos ont ot. (Jia7) 
In a similar manner, before low pass filtering, y,(t) is 
veau) eee na anf te (oar) 
After low pass filtering, y(t) would be known as the 
in phase component of the baseband noise, and y,(t) 
would be known as the quadrature component. 

At this point, we note that we have two baseband 
processes which have been generated from the bandpass 
received signal process, r(t). These two processes are 

ro(t) = s(t) + yo(t) = a/2 + y(t) (3.9) 
and 


ro(t) = s(t) + y(t) = y,(t). (3.10) 
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S2nee ty) contains ne information about the transmitted 
signal, we would like to disregard it and perform a 
standard baseband correlation receiver operation upon 
met ) only. It can be shown? that if the spectrum of 
y(t) is even over the bandwidth, then y,(t1) and Yalto) 
are uncorrelated for all observation times t, and to. 
Furtner, if y(t) were a gaussian process, then y,(t,) 

and yo(to) would be independent. Then r.(t) would 
contain no additional information about either the signal 
or the noise, and it could be disregarded with no loss 
meoptimality. 

Bandlimited low frequency atmospheric noise is 
definitely non-gaussian, but Feldman! hes noted that it 
has a spectrum which is determined solely by the band- 
limiting filter's frequency response squared, which we 
assume to be even over the bandwidth. Therefore, we can 
state that ¥_(ty) and Yglto) are indeed uncorrelated for 
all time instants t4 and to. Moreover, we assume that 
the statistical dependence between Yo(t,) and y (to) is 
sufficiently small so that we can ignore r,(t). We 
feel that this is not a limiting assumption, especially 
when viewed in the light that it allows us to simulate 
a baseband receiver as opposed to the nearly impossible 
task of simulating the passband receiver. We also invoke 


the fact that we are comparing the results of two 
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receivers under the same set of conditions to justify our 
assumption that ignoring r.,(t) is of little consequence. 

Now that we have transformed our received signal into 
a baseband signal, we will proceed with the sinulation 
of a baseband receiver. The received signal is now 
represented by equation 3.9. As a practical matter, we 
generated y,(t) on the computer by sampling y(t) at twice 
its bandwidth and multiplying the samples by the value 
of cos ernf,.t at the sample times. In the next section, 
we will state how we tested our receiver and noise model 
implementations on the computer by calculating the 


probability of error in the gaussian noise case. 


~2 Receiver Imnlementations 





In the last section, we described the menner in 
which we converted an untractable passband receiver 
Simulation problem into a more tractable baseband 
Simulation. In this section, we will present the digital 
model for our receivers. We will verify our single 
channel receiver by testing it with gaussian noise, 

We will check the simulation results with analytically 
calculated results, and we will show that these results 
verify our digital computer implementation of the 
Single channel receiver as correct. We will then test 
our model for the generation of the atmospheric noise 


Samples by comparing tne results in the nearly gaussian 
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quiet noise case with the gaussian results. These results 
will be similar if our implementation of the noise 
generator is correct. Since the noise generator is the 
same for all noise conditions (quiet, quiet night, tropical, 
and frontal), with the different conditions simulated by 
changing parameters only, we will then assume our noise 
generator is correct for all noise conditions. 

The basic single channel correlation receiver 
model is shown in computer flow graph form in Figure 3.1. 
We note that we have used a rather simple integration 
algorithm, which is a slightly modified trapezoidal 
rule. We ee as it as mee: 


1 i x(t) at x1 = (x(n) + x(n-1)) (3.11) 


en 
2 


where At is the interval between samples, NAt = T, and 
x(n) is a donee of x(t) at time n/t. We note that, in 
our case, we are integrating over one bit length, and 

we can only include N samples in our summation. This 

is because the sample called x(N) for the i" vit is 
also the sample called x(0) for the ae wat) Since 
the signal on each of these bits may be different, we 
must assign the end voint to one bit only. We therefore 
consider only samples x(0) through x(N-1) for each bit, 
and, in order to simplify our integrator, we assume that 
x(N) would aoe equal x(0). Thus, 


N- 
al Se Coeteaci jest = x(n), (3.12) 
T n=1 2 ~ N n=0 
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Figure 3.1 Single Channel Receiver 
Computer Flow Graph 
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which makes our integrator 
iD 
a f° x(t) at 1 
Tt 36 N 


which is the slightly modified trapezoidal integration 


Ned 
= x(n), (Gs ake) 
no 


rule that we have used. We feel that when N is relatively 
large (for example, 40, as in our case), the errors 
introduced by this simplification are small. 

In Figure 3.2; we show the results of a simulation 
of our single channel in which the noise was a gaussian 
process. We conducted three simulation runs of 1000 bits, 
at each of two bit lengths, averaged the probability of 
error estimates for each bit length, and plotted them 
versus the signal-to-noise ratio (SNR). We also 
calculated the probability of error for each SNR simulated, 
and those curves are plotted along with the corresponding 
Simulation results. AS can be seen, the calculated and 
experimental curves are in close agreement. We therefore 
conclude that our digital computer implementation of the 
Single channel receiver is correct. 

In Figure 4.4, we present the results of our 
Simulation runs using our atmospheric noise generator with 
the quiet noise parameters. We note that, as we expected, 
the results in this case are nearly the same as in the 
gaussian noise case, for the single channel receiver. 


We feel that these results verify our implementation of 
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aay 

the atmospheric noise model in the quiet noise case, and 
hence in all other cases, since they involve simple 
parameter changes in the generator. 

Finally, Figure 3.3 shows a computer flow graph of 
the dual channel receiver. We see that it is basically 
the same as the single channel implementation, with the 
exception that it incorporates the welehting CUNnCELOn., 
which is based upon the second channel noise envelope 
samples, inside the integration routine. We feel that, 
Since our results, shown in Figure 4.1 through Figure 
4.4, indicate significant improvements in performance in 
the more non-gaussian noise cases and comparable performance 
in the nearly gaussian quiet noise case, when compared 
with the single channel receiver, our dual channel 


receiver implementation is correct. 


3.3 Analog to Digital Calculations 


In testing our single chennel receiver versus 
analytically calculable results, we encountered the 
problem of relating eve enetion and digital noise power 
Spectra. We will calculate that relationship in this 
section. First, however, we will show how we calculated 


the probability of error per binary bit of information 


when the noise is Sosa eree 
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Define the average power in the information part of 
the received signal, r(t), as follows: 
Ba = Bey 2, (3.14) 
where T is the bit length and 


T 
Es =f s“(t) dt, (3s) 


where s(t) is the baseband signal portion of r(t), i.e. 
s(t) = + A/2. Define the source rate, R, as the rate 
at which binary bits are generated by the transmitter. 
Then, the energy per bit, Ey: is given by 

Ep = P /R. 


Therefore, 
T 2 
aay s“(t) at = (A/2)~. (3.16) 
pees (s/2) 
R 


The baseband SNR of the simulated signal and noise 
is defined as follows: 

SNR = E,/Nos (3517) 
where N)/2 is the variance of each independent sample of 
the gaussian noise process in the digital comvuter 
Simulation. However, in order to analytically compute 
the probability of error, we must obtain an expression 
for the power spectrum of a continuous white gaussian 
noise process, which is equivalent to the sampled noise 
process used in the simulation. We can do this in the 


following manner, in which n(.) represents the noise 





~ yy 
oe ess. 


Let the , bh noise sample be defined as 


AT 
(k) = it (t) dt, CAy160 
n 1S n 
then 
At - Gb 
B(n2(k)) 4 No/2 = BL J J, n(t,) n(t,) dt, dty) 
At At. 
= 1 ff B(n(ty) n(tg)) aty dts. (3.19) 
Ave Ona 
Recognizing that 
B(n(t,) n(ty)) = Neg §(ty-ta), (3.20) 
=e 


where §(.) is the Dirac delta function, we get 





No/2 = N,4/2s aazis) 
At 
and therefore, 
Neq = NoAt; (Be22) 
where At is the time interval between samples and Neq is 


the magnitude of the noise power spectrum of a continuous 
white gaussian process, analogous to the simulated process. 
Thus, the analog SNR corresponding to the digital SNR 
(equation 3.17) is given by 
SNR = EL/Neq = Ep/Nodte (3.23) 

For antipodal signalling schemes, the probability 
of error per bit, with gaussian noise, can be calculated 
as PO leneee. 


PCE) = Q( )/2E,/Neg’)s (3.24) 


where Q(x) is defined as 





‘ ~!}5= 
©0 z 
a(x) = f[ 1. e% / ax, (3.25) 


x 29 
whose complement, 1 - Q(x), is tabulated in various 
references (e.g. 8). Thus, in terms of the parameters 
of our digital simulation, the probability of error per 
binary bit, with white gaussian noise, can be analytically 


calculated by means of the Q(.) function as 


P(E) = Q( A 2E,/NoAt ee (2226) 


3.4 Chanter Summary 


In this chapter, we have explained why and how we 
converted our bandlimited receiver from a passband to a 
baseband system for digital simulation purposes. We went 
on to explain how we checked our implementations of the 
receiver and noise models by comparing simulated and 
analytical results using gaussian noise corruption and 
corruption by the quiet noise case of atmospheric noise. 
Finally, we gave the relationship between the simulated 
noise power spectrum and the analogous continuous power 
spectrum, which was necessary to calculate the probability 
of error per binary bit analytically. 

In the next chapter, we will present the results of 
our simulation runs using the atmospheric noise model 
described in Chapter 1. We will show the averaged 
results for flve runs of each of the four basic 
weather/noise conditions: quiet, quiet night, tropical, 


and frontal. We will show how we obtained our estimate 





Boe 
of P(E) per bit, and we will give bounds on that estimate. 
We will also give an example of the problems we 
encountered as we proceeded with the simulations. Finally, 
we will summarize the conditions under which our. results 
were obtained so that the reader may evaluate their 


usefulness, 
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Chapter 4 


SIMULATION RESULTS 


Before presenting our results, we will briefly 
Summarize our work up to this point. We discussed, in 
Chapter 1, bandlimited low frequency atmospheric noise, 
its causes, its non-gaussian nature, and a model for 


1 is 


this noise, The model, develoned by Feldman, 
especially suited for simulation on a digital computer, 
In Chapter 2, we discussed several receiver schemes, 
which could be used to decide which of two antipodal 
Signals, corrupted by atmospheric noise, was transmitted. 
Specifically, we discussed a linear correlation receiver, 
which is optimum when the noise encountered is gaussian, 
We then discussed a design, based upon the observed 
properties of the non-gaussian noise, which is considerably 
more complicated than the simple correlation receiver, 
but which should yield optimum performance in the 
presence of this noise. Finally in Chapter 3, we 
discussed the problems involved with a simulation of 
these receivers and bandlimited low frequency atmospheric 
noise, 

In this chapter, we present the results we obtained 
in our simulation of these receivers. We tested each 
receiver for each of four weather/noise conditions, 


which our noise model is capable of simulating, ranging 
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from the nearly gaussian quiet noise case to the extremely 
non-gaussian frontal noise case. YFror each receiver and 
weather/noise condition, we made five simulation runs, 
each run simulating the reception of 1000 bits cf binary 
informetion by the receiver. In order to characterize 
receiver performence, we formed an estimate of the 
probability of error per bit of received information for 
eacn rune We then averased the results of the five runs 
for each case and plotted these averages as functions of 
the signal-to-noise ratio of the received besebaend 
Signals. 

We present a discussion of the method in which we 
formed our estimates of probability of error, and we 
give bounds upon our estimating procedure. We also 
bound our results by presenting tables of certein 
Statistics of the results of the five individual runs 
for each case. We discuss briefly one problem we 
encountered during the simulation runs, and why this 
problem existed and how we solved it. Finally, in the 


last section, we present a summary of this research. 


4.1 Results and Bounds 

In Figures 4.1 thrceugch 4.4, we have plotted the 
results of our simulation runs for frontal, tropical, 
quiet night, and quiet weather/noise conditions, 


respectively. The curves shown are plotted on semi-log 
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Figure 4.1 Frontal Noise 
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Figure 4.2 Tropical Noise 
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Figure 4.3 Quiet Night Noise 
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paper, with the ordinate representing the probability of 
error per bit and the abcissa being the signal-to-noise 
ratio in decibels (20 log SNR). There are three curves 
shown in each figure. These curves represent the results 
of the single channel correlation reciever, the dual 
channel receiver, and the ad hoc receiver (described in 
Section 2.3) wnich we used as a check upon our dual 
channel receiver. 

The curves shown in each graph represent the average 
of five simulation runs of 1000 bits of received infor- 
mation each. The bit length for all runs was 20 msec, 
making the simulated lergsth of each run 20 seconds. The 
eel <tolnoise ratio scale was defined in the following 
manner; It is necessary to realize that the noise power 
level is vastly different from one weather/noise condition 
to any other. However, the pawn imeed low frequency 
atmospheric noise model used in these simulations consists 
of two linearly combined functions, a background low level 
gaussian process, and a process which represents the high 
energy burst behavior of the nodcene The background 
gaussian noise is assumed to have been generated in the 
receiver itself (receiver front-end oe and therefore, 
it is independent of the particular weather/noise condition 
being simulated. Since we would like our SNR scale to 


be independent of the weather/noise condition being 
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simulated, we base it upon the signal power and the power 
level of the background gaussian process. 

Thus, our SNR is given by equation 3.23, which is 
reveated here: 

SNR = E,/Nodts 
where Ey is the energy per bit, as given by equations 3.14 
through 3.16, No/2 is the variance of each independent 
sample of the background gaussien noise process of our 
noise model (all the parameters of our noise model have 
been normalized such that No/2 a ee and At is the 
Sampling interval. 

The estimate of the probability of error per bit 
was formed by counting the number of incorrect decisions 
made by the receiver and dividing that number by the 
totel number of bits simulated. For example, in a run 
where we simulated 1000 bits, and where the receiver 
made 351 incorrect decisions (we have assumed that the 
probability of making an error is independent of the 
Signal, either a binary one or a binary zero, which was 
sent), our estimate of the probability of error per bit 
in that run would be 0.351. 

We now present an evaluation of our estimation 
procedure, First, we assume that the probability that 
an error is made on any particular bit, for any 


particular receiver and weather/noise condition, is 
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constant end can be denoted by the value p. Then the 
probability that a correct decision is made on any bit 
is simply (1 - p). We assume that the decision made on 
each bit is independent of any, and all, other bits. 
Now we can assign a random variable, say Ks to the 
th 


outcome of the i aecision of the receiver. The random 


variable, k is defined to have the following probability 


4? 
density function: 


de) 0) 


i al Ds; (4.1a) 


P(k, = 1) =p. (4.10) 


Therefore, we see that k 


error is made on the ae bi Cgaend Ky is zero otherwise. 


i 1s one if, and only if, an 


Now we define the random variable s, such that 
N | 
s=1/N > ky, (4.2) 
ta 


where N is the number of bits simulated. We know that 
the set of Ky is a set of independent identically 
distributed random variables, each having the probability 
density function given in equations 4.1 above, which we 
note is the well defined binomial density. Therefore, 
the mean and variance of each ky are given Belowi” 

E(k,) =p, (4.3) 

Var(k,) = p (1 - p) 

=p- p*. (4.4) 


Now, we see that, by the central limit theorem, if N is 
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large, the probability density function for s is gaussian 


with mean and variance described by 


N 
E(s) = E(1/N = k,) 
i=l 
N 
= 1/N 5 E(x,) 
i=1 
=a (4.5) 
and 
Var(s) = E(s“) - (E(s))* 
NN 
= E(1/Né = = Ky ky) = pé 
i=1 j=1 
= (p - p*)/N, ee) 


where we have used the independence of the Ky We make 
note of the facts that s is an unbiased estimate of p, 
and that the variance of s decreases toward zero as the 
number of bits simulated increases... 

Finally, in order to establish a bound upon our 
estimates of probability of error per bit, as given by 
Ss in equation 4.2, we used a reo of the complement 
of the Q(.) function, described in equation 3.25. We 
note from the table that, with 95% confidence, our 
estimate of the probability of error per bit will be 
within two standard deviations of the actual value. We 
also note that the variance of our estimate, and hence 
its standard deviation, is greatest when the actual 


value of p is 0.5. We therefore use this case for a 
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worst case analysis of our estimation procedure. 
When p = 0.5, the variance of s is 
Var(s) = (0.5 )-00"e5)/N =f02-> |. (4.7) 
In our case, we have made five runs of 1000 bits each; 
thus meking N = 5000. Therefore, 
Vamsi = (0mes7 5009R (4.8) 
and the standard deviation, @,, is 
a = ,{ 0.25/5000' = 0.00708, (4.9) 
and thus, 
20, = 0,01416,. (4.10) 
From the above, we can see that, with 95% confidence, 
our estimates of the actual probability of error per 
binary bit should be within less than + 0.014 of the 
actual values. As a matter of practical concern, we note 
that each run of 1000 simulated bits of information was 
conducted independently of the others. We therefore 
also present a summary of the results of the five runs 
for each weather/noise condition and receiver. These 
results are presented in tabular form in Tables 4.1 
through 4.4. In these tables, we give the average of 
probabilities of error of the five runs for each case and 


the maximum deviation from that average. 


Hee Diffileniey 


represen cence es 


The most serious difficulty encountered in this 


research was an apparent breakdown in our single channel 
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Table 4.4 Quiet Noise 
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linear correlation receiver model. Our original intention 
was to make only one simulation run, of 1000 bits, for 
each receiver and weather/noise condition while varying 
both signal magnitude and bit length. We note that an 
increase in either signal magnitude or bit length results 
in an increase in SNR, and therefore should result in a 
decrease in the probability of error per bit. For our 
dual channel receiver, this was indeed the case, but with 
the single channel model, we found that the probability 
of error was increasing slightly with increasing bit 
length for all signal magnitudes below a certain, non- 
constant threshold. 

We reasoned that it would be interesting to see how 
well each receiver performed as a function of both the 
Signal magnitude and the bit length. Therefore, we 
thought that we were keeping all eeen factors constant 
and varying only the factor of interest, either the bit 
length or the signal magnitude, from run to run. We 
must diverge here slightly and exnvlain the manner in 
which we generate the noise process upon a digital 
computer. To generate any random process upon a digital 
computer, we in fact are generating a sequence of 
pseudorandom numbers, these numbers being samples of the 
noise process. In order to generate this sequence, we 
must initialize all random number generators in the 


computer program (see Figure 1.1). Furthermore, if we 
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initialize all the random number generators the same way 
in every run of the same program, we will renerate the 
identical random number sequence.” 

Now we can see that when we varied the signal 
magnitude only, and kept the bit length constant, the 
noise process which corrupted each bit was the same for 
all signal magnitudes since we initialized all random 
numbers in the same manner at the start of each run. 
However, when we varied the bit length, while keeping 
the signal magnitude constant, the noise process seen by 
each bit changed radically with bit length, regardless 
of how we initialized the random number generators. Thus, 
we were varying two factors, both the noise process and 
the bit length. We made several more runs, with different 
initial conditions upon the random number generators, and 
we found that the probability of error did indeed decrease 
as we increased bit length in these runs. We therefore 
reasoned that our first run had been made with a none- 
typical noise process that caused the probability of 
error to increase slightly with increasing bit length. 

We subsequently decided that we should make several 
runs, starting with different initial conditions upon 
the random number generators each time, for each bit 
length and each signal magnitude. However, due to time 
constraints, we found that we could only make several 


runs for each weather/noise condition, receiver, and 
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Signal magnitude at one bit length. We feel that in this 
way our results are typical of the actual results to be 
expected for that bit length. We chose a bit length of 
20 msec, about the shortest bit length practical in low 
frequency digital communication systems, 2 because it 


resulted in the shortest time per simulation run. 


4.3 Summarv of Reseerch Conditions 





In this section, we present a summary of the 
receiver and noise models we used in the simulation 
runs which we made. We also give a short summery of 
the results we obtained. 

The model for bandlimited low frequency atmospheric 
noise which we used was developed by einen. 2 This noise 
model accurately models the first order statistics and 
the time structure of the actual bandlimited noise 
waveform for four weather/noise conditions. These four 
conditions include the nearly gaussian quiet noise case, 
the quiet night noise case, the tropical noise case, and 
the most severe frontal noise case. These four conditions 
model noise which is typical of large of areas of the 
world over significantly large periods of time. The 
quiet noise conditions exist when the atmsophere in the 
vicinity of the receiver consists of a stable alr mass. 
The tropical conditions exist when the air is wam, moist, 


and unstable. Quiet night conditions represent a 
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transitional period between quiet and tropical conditions. 
And, of course, frontal noise conditions, the most violent 
conditions, exist when there is a front, with much 
thunderstorm activity, in the vicinity of the receiver. 

A canonic noise generator, useful in comvuter simulations, 
and a table of associated parameters are given in Figure 
1,1 and Table 1.1 of this work. 

For our simulation, we chose binary antipodal 
Signalling, with a coherent PSK modem. The carrier 
frequency of our signals was 65 kHz and the receiver's 
bandlimiting filters had a bandwidth of 1 kHz. We also 
measured the noise envelope on a second frequency channel, 
with bandwidth 1 kHz, centered at 83 kHz. These frequencies 
were chosen because they are typical of digital communication 
Systems, and we have the parameters for simulation of 
the noise model at these frequencies.~ 

As a standard of comparison, we simulated a single 
channel linear correlation receiver, using a 20 msec bit 
length, operating in the above described low frequency 
atmospheric noise. Our test statistic, an estimate of 
the receiver's probability of error per binary bit, was 
calculated simply by counting the number of incorrect 
decisions made by the receiver and dividing by the 
number of bits simulated. We made five runs, of 1000 
bits each, for each weather/noise condition for various 


signal-to-noise ratios. 
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We then tested, in a similar manner, a dual channel 
receiver which uses the fact that, above a certein 
threshold, the value of the noise envelove on the 83 kHz 
channel and the value of the noise waveform on the signal 
channel are statistically dependent. Due to the sampling 
rates involved, and limitations of the noise model, we 
simulated both of the above receivers in the baseband 
rather than the passband. 

We found that this dual channel receiver performed 
better than the single channel receiver for the three 
worst noise conditions (frontal, tropical and quiet night), 
and as the signal-to-noise ratio increased, the performance 
improvement exceeded an order of magnitude. However, in 
the quiet noise case, we found that the dual channel 
receiver yielded approximately the ane performance as 


the single channel model. 


4.4 Chanter Summary 


In this chapter, we have presented our results of 
the simulations of the single and dual channel receivers. 
We also presented an analysis of the estimate of probability 
of error per bit which we used, and we gave bounds upon 
our results. We explained a difficulty we encountered 
in the course of the simulations, and we explained how 
we overcame it. Finally, we presented a summary or the 
research reported in this work. In the next chapter, we 


present our conclusions and ideas for further research. 
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Chaepver: 5 


CONCLUSIONS AND SUGGESTIONS FOR FURTHER RESEARCH 


We first present a brief summary of this research. 
In Chapter 1, we described low frequency atmospheric 
noise and a model for this noise. In Chapter 2, we 
presented two receiver designs: a linear correlation 
receiver, the optimum receiver when the noise corrupting 
the signal is gaussian; a dual channel design, based 
upon the ooserved non-gaussian characteristics of 
bandlimited low frequency atmospheric noise, which 
should be optimum for tnis noise. In Chapter 3, we 
presented the calculations necessary to relate the 
digital simulation to the analog system being modelled. 
In Chapter 4, we presented our results for the two 
receivers described above, and some bounds upon those 
mesults. In this chapter, we will give our conclusions 
based upon our results, and we will also give some 
suggestions for further research with the dual chan- 


nel receiver design, 


5-1 Conclusions 


Our results, shown in graphical form in Figures 
+e1 through 4.4, were favorable. We found that the dual 
channel receiver yielded significant performance improve- 


ments over the single channel receiver in the three more 
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violent weather/noise conditions (frontal, tropical, and 
quiet night). In the nearly gaussian quiet noise condition, 
the dual channel receiver performed avnproximately the 

same as the single channel linear receiver. The 
performance improvements in the more violent noise cases 
were more than an order of magnitude as the signal-to- 

noise ratio increased. 

Tne results were as we expected. We note that the 
Sinsle channel linear correlation receiver is the optimum 
gaussian receiver, and it tends to be optimum in the 
quiet noise case. However, as the noise becomes more 
non-gaussian, we note that the dual channel receiver 
performs much better. Finally, in the frontal noise 
case, the most violent noise case, the single channel 
receiver could do little better tnan guess, while the 
duel channel receiver made correct decisions with about 
98% confidence at the same signal-to-noise ratio. 

We therefore conclude that the dual channel receiver, 
while beins approximately twice as complicated as the 
Single cnannel receiver, yielded a significant improve- 
ment anid should be investigated further. It must be 
noted tnat we have simulated only the simplest case, 
using coherent PSK signals, but we feel thet our results 
indicate that the dual channel receiver design is a 
potentially powerful design for low frequency digital 


communication systems. 
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e2 Suerestions for Further Research 





We have tested our dual channel receiver design 
against the optimum receiver for gaussian noise, but 
we have not tested it against a receiver which is more 
similer to those currently used to receive signals 
corrupted by non-gaussian atmospheric noise. Practical 
correlation receiver designs for use in low freauency — 
signalling systems use nonlinear devices, placed before 
the integrator, to limit the effects of the non-gaussian 
high energy noise bursts. Typical nonlinear devices 
which ere used, or suggested for use, are clippers, 
limiters, and hole Puncher aaa (see Figure 5.1) 

We suggest thata possible research program, based 
upon our results, would be to determine the optimum 
nonlinear device to be placed in ee single channel 
receiver, and to determine the optimum point at which 
limiting, or attenuating, occurs. Simulations could be 
run, using this ootimum nonlinear receiver, and the 
results compared with the results we have obtained with 
the dual channel receiver. Since the single channel, 
nonlinear receiver would be simpler to implement than 
the dual channel design, it would be interesting to 
determine how much performance improvement, if any, is 
possible with the dual channel design over the single 
channel nonlinear desien. 


In our research, we assumed we had phase coherent 
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Figure 5.1 Nonlinear Devices 
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receivers, and we postulated tnat the results we obtained 
would give an indication of the performance to be expected 
in the case of incoherent receivers. Several standard 
communication theory references (e.g. 3,+) could be used 
to find inplementations of incoherent receivers, and the 
test we made could be redone using these incoherent 
receivers. 

We found several references to an incoherent receiver 
en which uses a combination of energy detection 
anc 2 nilot tone to determine the nhase of the received 
Sicnal. This receiver has been simulated and results 
have been presented for the case in which the corrupting 
noise is eaussian.” It would be interesting to imvolement 
this receiver in a dual channel design, and to test it 
in the vresence of the non=-gaussian atmospheric noise. 

The information from the second channel can also be 
used in an error detection scheme. Since the noise has 
a time varying power level, the actual signal-to-noise 
ratio for each bit of received signal may rluctuate. The 
second channel noise envelope measurements can be used 
to determine the actual signal-to-noise ratio on any bit, 
and when tne ectual SNR is below a certain threshold, 
the receiver cen indicate that the bit is probably in 
erro 

Since we found that the single cnannel receiver 


is indeed the optimum receiver during the quiet noise 
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periods, while the dual channel receiver performs much 
better in the more severe noise periods, it would be 
convenient to have some device which determines when 
to switch between the single end dual channel receivers. 
Furthermore, the thresholds for the dual channel design 
very with the particular weather/noise condition in 
existence. Therefore, it would be beneficial to Gesign 
an algorithm which samples the noise envelope on the 
second channel and determines which receiver and/or 


which threshold (v,5 from equation 2.16) value to use. 
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